Finite-temperature properties of hard-core bosons confined on 
one-dimensional optical lattices 
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We present an exact study of the finite-temperature properties of hard-core bosons (HCB's) 
confined on one-dimensional optical lattices. Our solution of the HCB problem is based on the 
Jordan- Wigner transformation and properties of Slater determinants. We analyze the effects of the 
temperature on the behavior of the one-particle correlations, the momentum distribution function, 
and the lowest natural orbitals. In addition, we compare results obtained using the grand-canonical 
and canonical descriptions for systems like the ones recently achieved experimentally. We show that 
even for such small systems, as small as 10 HCB's in 50 lattice sites, there are only minor differences 
between the energies and momentum distributions obtained within both ensembles. 

PACS numbers: 03.75.Hh, 05.30.Jp 



I. INTRODUCTION 

The study of ultracold quantum gases loaded on opti- 
cal lattices has become a very active area of experimental 
and theoretical reseach in recent years. Optical lattices 
enable enhancing interactions between atoms in weakly 
interacting Bose-Einstein condensates (BEC's) and re- 
ducing the effective dimensionality of the system. They 
allow the experimental realization of strongly correlated 
bosons, well described by the Bosc-Hubbard Hamiltonian 
HI j with the consequent observation of the superfluid- 
Mott-insulator transition 0, Q ■ In addition, optical lat- 
tices have been used to obtain one-dimensional (ID) sys- 
tems @ and to examine the superfluid-Mott-insulator 
transition in ID 

Due to the strong effects of quantum fluctuations 
and the possibility of obtaining exact theoretical results, 
ID systems are a very attractive laboratory for both 
experiments and theory. Theoretically, it was shown 
by Olshanii that in ID in regimes of large scattering 
length, low densities, and low temperatures bosons be- 
have as a gas of impenetrable particles known as hard- 
core bosons (HCB's) @. Such a ID gas (also recently 
called a Tonks-Girardeau gas) was introduced by Gi- 
rardeau, who established an exact mapping between 
these strongly correlated bosons and noninteracting spin- 
less fermions ^(J- Since then ID HCB's have been ex- 
tensively studied b y d ifferent techniques in both homo- 
geneous: Jnl EH fl3L ll4L ITU Hr3 | and harmonically trapped 
[ill ITiL Il9f systems. 

The experimental realization of ID HCB's followed af- 
ter more than 40 y ears of the theoretical introduction of 

the model [UEII, with and without E3 an addi- 
tional lattice along the ID axis. The additional ID lat- 
tice [2(| facilitates the achievement of the HCB regime 
with respect to the continuum case. It allows experi- 
mentalists to change the effective mass of the particles 
and, consequently, the ratio between interaction and ki- 
netic energies [2(| . Although at very low densities (when 
intcrparticle distances are much larger than the lattice 
spacing) HCB's on a lattice are equivalent to HCB's in 



continuous space, this is not the case for arbitrary fill- 
ings . On ID lattices the HCB Hamiltonian can be 
ma ppe d onto the ID XY model of Licb, Schulz, and Mat- 
tis J23j. For periodic systems this model has been also 
studied extensively in the literature [24l I2M |26| . More 
recently, renewed interest has arisen on the properties 
of HCB's when additional confining potentials are intro- 
duced, as the case relevant to experiments |20l |. 

Remarkably, even in trapped inhomogeneous systems 
power-law behavior known from the periodic case is 
present (2l|. The one-particle density matrix exhibits a 
universal power-law decay with exponent —1/2 indepen- 
dent of the power of the confining potential [2l} . These 
quasi-long-range one-particle correlations generate qua- 
sicondensates with occupations scaling proportional to 
^/Nb (with Nb the number of HCB's in the system) [27j . 
The nonequilibrium dynamics of HCB's on ID lattices 
has also been shown to display very interesting features. 
Quasicondensates of HCB's emerge at finite momentum 
when the system starts its free evolution from a pure 
Mott-insulating (Fock) state [2^|. In adition, it was 
shown in Rcf. [2!j that in ID when there is no Mott 
insulator in the trap, the momentum distribution of ex- 
panding HCB's rapidly approaches that of noninteracting 
fermions p9l |. 

In this work we present an exact study of the finite- 
temperature properties of HCB's confined on ID optical 
lattices. Following the spirit of Refs. [2H l2Sj . we de- 
velop an exact numerical approach based on the Jordan- 
Wigncr transformation, which maps HCB's on a lattice 
onto noninteracting spinless fermions. We will focus on 
the effect of temperature on the off-diagonal behavior of 
one-particle correlations and related quantities like the 
momentum distribution function and the natural or- 
bital occupations. The natural orbitals (<p v ) arc defined 
as the eigenfunctions of the one-particle density matrix 



A' 



(1) 



3=1 

and have occupations X^. (They resemble one-particle 



2 



states in these strongly interacting systems.) In dilute 
higher-dimensional gases, when only the lowest natural 
orbital (the highest occupied one) scales ~ Nb, it can be 
regarded as the BEC order parameter |3l| ■ Here we will 
show that in ID even at very low temperatures (T) , when 
the energy of the system is almost identical to the one at 
T = 0, the momentum distribution and the lowest nat- 
ural orbital occupations can exhibit significant changes 
with respect to their values in the ground state. 

The exposition is organized as follows. In Sec. II we 
describe our exact approach to study finite-temperature 
systems. In Sec. Ill we discuss the properties of HCB's in 
a perfect box (an open system). HCB's confined in har- 
monic traps are analyzed in Sec. IV. Since in this work 
we follow a grand-canonical approach to study finite- 
temperature properties, in Sec. V we compare exact re- 
sults obtained from a grand-canonical calculation with 
results obtained from a canonical one for small lattice 
sizes, like the ones recently achieved experimentally |2f| . 
Finally, the conclusions are presented in Sec. VI. 



II. EXACT FINITE-TEMPERATURE 
APPROACH 

In this section we detail the exact approach followed 
to study the finite-temperature properties of HCB's con- 
fined on ID lattices. The HCB Hamiltonian can be writ- 
ten as 



i7 = -*E(XVi + H 



Vi^a^m, (2) 



with the additional on-site constraints 



6? 







(3) 



which avoid double or higher occupancy. The bosonic 
creation and annihilation operators at site i are denoted 
by b\ and b t , respectively, and the local density operator 
by n, = b\b i . The brackets in Eq. J2J) apply only to on- 
site anticommutation relations; for i ^ j, these operators 
commute as usual for bosons = 0. In Eq. J5J), 

the hopping parameter is denoted by t and the last term 
represents a harmonic trap with curvature V2. 

In order to exactly calculate HCB properties, we use 
the Jordan- Wigner transformation [32] 



/3=1 



t — ± 

k=l[e i * f ^f i , (4) 



which maps the HCB Hamiltonian onto the one of non- 
interacting spinless fermions, 



x? n{, (5) 



where fj and /j are the creation and annihilation oper- 

ators for spinless fermions at site i and = fl f i is the 
local particle number operator. 



The mapping as presented above is only valid for open 
systems, as relevant for confined bosons in experiments 
[2], il|. In such cases HCB's and fermions have ex- 
actly the same spectrum. In order to deal with ID cyclic 
chains, with N lattice sites, one needs to consider that 



N 



b l b N = -/l/jV ex P \ i7T Yl 



(6) 



0=1 



so that when the number of particles in the system 

[J2i( n i) — J2i( n i) = Nb] is odd, the equivalent fermionic 
Hamiltonian satisfies periodic boundary conditions; oth- 
erwise, if Nf, is even, antiperiodic boundary conditions 
are required in Eq. (J5J. 

Since for finite temperatures we will consider a grand- 
canonical ensemble — i.e., a system with fluctuating num- 
ber of particles — in order to avoid the dependence of the 
equivalent fermionic Hamiltonian on Nb we restrict our 
analysis to the open case. In this case the nontrivial dif- 
ferences between the properties of HCB's and fermions 
are only in off-diagonal correlation functions. 

For finite temperatures, and within the grand-canonical 
formalism, the HCB one-particle density matrix can be 
written in terms of the equivalent fermionic system as 



Pir- 



— Tr < b\b, exp 



H - n l I A'bT 



1 f J 1 

[ 0=1 



(7) 



x exp 



H F — /j, n{ J /k B T 



^ cxp(— iirnf) 

7=1 



where, in addition to Eqs. Q and JSJ, we have used the 
cyclic property of the trace. In Eq. fi denotes the 
chemical potential, k B the Boltzmann constant, T the 
temperature of the system, and Z the partition function 



Z = Tr- 



j exp - ^H F -^J2 n ij l kBT } 



(8) 



To calculate traces over the Fock space we will take 
advantage of the fact that in the equivalent fermionic 
system Fock states are Slater determinants, 



N f N 



1=1 i=i 



(9) 



with Nf the number of fermions and 



/ Pu P12 
P21 P22 



\Pm Pn2 



PlN f \ 



P 



2N f 



(10) 



PNN f / 
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the matrix of the components. 

The action of exponentials bilinear on fermionic cre- 
ation and annihilation operators, as the ones on Eqs. (0 
and JSJ, on Slater determinants generates new Slater de- 
terminants [33|,|3 

/ \ N f N 

exp £ fjXijfj = II E Pl 3if) 1°)' ( u ) 



where 



P' 



(12) 



Using this property one can prove the following iden- 
tity for the trace over the fermionic Fock space |U E3 



Tr 



ij ) \ kl ) 

exp I fL Z mnf n j 



det [I + 



e x e Y 



(13) 



which immediately allows one to calculate the partition 
function as 



Z 



det [i + e -(H F -„D/fcBT 

JJ [l + e -(£«-M)ABT 



(14) 



where I is the identity matrix. The last equality was 
obtained after diagonalizing Hamiltonian ijjjj, HfU = 
UE, U is the orthogonal matrix of eigenvectors, and E 
is the diagonal matrix of eigenvalues. 

The trace in Eq. J7J) is calculated along the same line. 
For i ^ j , we notice that 



fifj = CX P ( E frn A rnnf n j ~ 1, 
\ mn / 



(15) 



where the only nonzero element of A is Ay = 1. Then, 
for i =/= j, pij can be obtained as 

Pij = - ( det f 1 + ( r + A)OiUe- (E ~ flI)/fcsT U t 2 



- det 



I + OiUe~ (B -' a )/* aT U t 02] } ■ (16) 



Oi (O2) is diagonal with the first j — 1 (i — 1) elements 
of the diagonal equal to —1 and the others equal to 1. 

The diagonal elements of the one-particle density ma- 
trix are the same of noninteracting fermions [see Eq. (J7J 
for i — j] and can be easily calculated as |33l. l34l| 



j + e -(H P -^i)/fc B T 



U I 



-(E-ftI)/fe B T 



ut 



(17) 



As usual, the chemical potential is fixed using the relation 
Nb = pa to obtain the desired number of particles in 
the system. 



III. HARD-CORE BOSONS IN A BOX 

In this section we study the finite-temperature prop- 
erties of HCB's on a perfect box. In this case the HCB 
Hamiltonian can be written as 



H = —t 



N-l 

E 



b]b i+1 



H.c. 



(18) 



with the additional on-site constraints (j2J). 

The above Hamiltonian is particle- hole symmetric, like 
the one of periodic systems, under the transformation 
hi = b\, h\ = hi {h\ and hi are hole creation and anni- 
hilation operators). The particle- hole symmetry implies 
that the off-diagonal elements of the one-particle density 
matrix for N b HCB's [pij(N b )] and for (N - N b ) HCB's 
[pij(N — Nb)] are identical. Diagonal elements satisfy 
the relation pa(Nb) = 1 — pu(N — Nb). This leads to a 
momentum distribution function 



1 E' 

N ^ 



-ik(xj —xi) 



Pjh 



which satisfies the relation 



n k {N b 



-k{N-N b )+[ 1 



N -N b 
N/2 



(19) 



(20) 



In contrast to periodic systems where the natural or- 
bitals [Eq. (J)] are momentum states 0, |2tJ , this is not 
the case in a box. (The system is not translationally 
invariant.) In Fig. ^ we show the lowest-natural-orbital 
wave function in a box at different temperatures. We 
have normalized it as 



ip = i? 1 / 2 0°, R = (NbN) 1 ^ 2 , 



(21) 



so that ip° vs x/N is independent of the system size when 
the density p = N b /N is kept constant. 




FIG. 1: Wave function of the lowest natural orbital at differ- 
ent temperatures in half-filled systems with 1000 lattice sites. 

Figure H shows that at finite temperatures the weight 
of the lowest natural orbital increases in the center of the 
system, departing from the constant value it would have 



4 



in the periodic case (k = state). Still, we find that 
qualitatively (and quantitatively) the natural orbital oc- 
cupations behave very similarly to the occupations of the 
momentum states so that for the box we will restrict our 
analysis to rife. The natural orbitals will be relevant to 
the discussion in the harmonic trap where their behavior 
can be qualitatively different to the one of rik- 

In Figs. Ot a )~Et c ) we show the HCB momentum dis- 
tribution function for half-filled systems with N = 1000 
and different temperatures. We have plotted as dashed 
lines the ground-state results for comparison. The ef- 
fects of small but finite temperatures are dramatic. This 
can be better seen in Figs.Efa) and[5{b) where the ener- 
gies of the finite-temperature systems are almost identical 
to the ones of the ground state. For UbT = O.Olt, the 
relative energy difference (5E = [E(T) - E(0)]/\E(0)\) 
between the finite-temperature system [E(T)] and the 
ground state [E(0)} is 6E ~ 0.01%. In Fig.^a) one can 
see that the k = momentum peak is already around 
2/3 of the one at zero temperature. For the case in Fig. 
EJb), 5E ~ 0.4% and the peak at n^, =0 has already re- 
duced almost 5 times. At fc^T = 0.5i in Fig. (a), the zero 
momentum peak has practically disappeared. 




FIG. 2: Momentum distribution function of HCB's (a)-(c) 
and noninteracting fermions (d)-(f) for half-filled systems 
with 1000 lattice sites at different temperatures (solid line). 
The temperature [energy] of the system in each case is ksT — 
O.Olt [E= -6.362 xl0 2 i] (a),(d), O.lOf [-6.336 x 10 2 t] (b),(e), 
and 0.50t [-5.585 x 10 2 t] (c),(f). The dashed line in all the 
figures depicts the ground-state result (E — —6.363 x 10 2 t). 

As opposed to the HCB momentum distribution func- 
tion, we have plotted in Figs. Eld) EI0 the momen- 
tum distribution function of the equivalent noninteract- 
ing fermions. These figures not only show the differ- 
ences between the shape of the momentum distributions 



in both cases, but also the fact that they are affected very 
differently by the temperature. In the fermionic case it is 
well known that the changes on rik occur only around the 
Fermi surface and are of order fcgT, so that in Fig. Eld) 
one cannot notice the differences between the finite- and 
zero-temperature cases. In Fig.EJe) they are very small, 
and only when fc^T becomes of the order of t [Fig. Elf )] 
can one see a large deviation of the finite-temperature nk 
with respect to the one in the ground state. 

The zero-temperature peaks in the HCB rik [Figs.Efa)- 
Figs. Et c )] reflect the presence of quasi-long-range one- 
particle correlations [24l l25l l2a l27j : i.e., there is a power- 
law decay pij ~ \xi — Xjl -1 ' 2 . In these ID systems any 
finite temperature generates an exponential decay of p^ , 
which destroys the quasi-long-range correlations present 
in the ground state. This exponential decay is the one 
producing dramatic effects in nk- 

In Fig. El we show the decay of one-particle correlations 
for the same systems of Fig. El At very low temperatures 
(ksT = O.Olt) the one-particle density matrix follows the 
ground-state result over a certain distance, which reduces 
with increasing the temperature, approximately up to the 
point where the exponential decay sets in. Our results in 
Fig. El can be compared with the ones obtained by other 
means for (SfSj) in the ID spin-1/2 isotropic XY model 
|35j . to which HCB's can be mapped. Apart from a (1/2) 
normalization factor the results agree. 




100 200 300 400 500 



FIG. 3: Decay of the one-particle density matrix at differ- 
ent temperatures in half-filled systems with 1000 lattice sites. 
Thin solid lines following finite-temperature results exhibit 
exponential decays, while the one following the T = result 
exhibits a power-law decay ~ xj 1 ^ 2 . In all cases we measured 
pij fixing j in the middle of the box — i.e., Xj = 

The quantity of relevance to characterize the finite- 
temperature exponential decay of the one-particle den- 



sity matrix ~ e 



(Fig. EJ is the correlation 



length £. This quantity is of experimental importance 
since for £ > iV the HCB gas (essentially) exhibits at fi- 
nite temperatures properties of the ground state. At low 
temperatures, fc^T < t, the correlation length decreases 
as £ ~ 1/T with increasing temperature. This is shown 
in Fig. El 



5 



1000 




100 



FIG. 4: Correlation length vs temperature (in units of t) in 
half-filled systems with 1000 lattice sites. We have plotted £ 
as a thick solid line and the second moment of the one-particle 
density matrix £ as a dashed line (see text). The thin solid 
line is the result of a fit f/a = l.27t/k B T. 



A way of seeing the effects that a finite-temperature 
correlation length produces in these bosonic systems is 
to study how the occupation of the zero-momentum state 
scales with the number of particles (or the system size) 
when the density is kept constant. Results for nk=o vs 
Nb are presented in Fig. |SJ There we have plotted re- 
sults for as many temperatures as in Fig. El so that one 
can see at what system size the finite-temperature re- 
sults depart from the ones of the ground state. Since the 
system size is twice the number of particles (they are at 
half filling), one can then notice, with the help of Fig - EI 
that the mentioned departure indeed occurs for system 
sizes larger than the correlation length. For example, a 
half- filled box with 20 HCB's (a filling similar to the one 
achieved experimentally in Rcf. 20]) would have a mo- 
mentum distribution function very similar to the one in 
the ground state up to a temperature fc^T = 0.05i. At 
zero temperatures rik=o scales proportionally to ^/IV^; i-e, 
it diverges when Nb — ► oo, reflecting the power-law decay 
of one-particle correlations shown in Fig. 13 23. 
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FIG. 5: Scaling of rife = o vs Nb for a constant density p = 0.5 
and the same values of the temperature chosen in Fig. [J] 



The one-particle correlation length not only depends 
strongly on the temperature, but also on the density in 
the system. (In Fig. El we have only shown results for 
the half-filled case.) The dependence of the correlation 
length on the density, for two values of the temperature, 
is depicted in Fig. El Notice that both curves are symmet- 
ric with respect to p = 0.5 due to particle-hole symmetry. 




FIG. 6: Correlation length vs density for systems with N = 
500, k B T = O.lf (solid line) and k B T = 0M (dash-dotted 
line). Thick lines depict £ and thin lines the second moment 
of the one-particle density matrix £. 

The strong dependence of the correlation length on the 
density represents a difficulty for defining £ in inhomo- 
geneous systems, like the ones achieved experimentally 
where HCB's are trapped in harmonic confining poten- 
tials. (In a box the density is not exactly constant, away 
from half and integer fillings, due to Friedel oscillations, 
but they reduce with increasing system size.) An alter- 
native definition to the correlation length £ may be given 
as the second moment of the one-particle density matrix 



\ 



1 Ey 0* 



v 2 

Xj ) pij 



Eij Pij 



(22) 



In Fig. El we have plotted £ along with £. When 
1 < £ -c N both £ and £ are very similar. For the 
lowest temperatures, in Fig. El we considered systems 
with 1000 lattice sites, which are not much larger than 
the correlation length. This is the origin of the differ- 
ences between £ and £ observed for large values of £. At 
high temperatures (ksT > t) the value of £ is completely 
dominated by the very-short-distance sector of the one- 
particle density matrix, so that £ and £ are expected to 
be very different. At intermediate temperatures one can 
use £ as a good estimate of £. 

In Fig. El we have also plotted £ along with £ so that 
one can realize how the inclusion of the short-range part 
of the one-particle density matrix in £ produces different 
effects for low densities, where £ > £, and high densities, 
where £ < £. Still the overall behavior of £ is similar to 
the one of £. In the next section we will rely on £ for 
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estimating the correlation length in harmonic traps and 
also for comparing it to the one in the box. 



IV. HARD-CORE BOSONS IN HARMONIC 
TRAPS 

We study in this section HCB's trapped in harmonic 
potentials. The addition of a confining potential gener- 
ates a position-dependent density profile where, at zero 
temperature, supcrfluid and Mott-insulating regions can 
coexist. In the next two subsections we analyze the ef- 
fects of the temperature on density and momentum pro- 
files of systems in which the ground state is (i) superfluid 
fSec. HV A"|) and (ii) a coexistence of supcrfluid and Mott- 

In Sec. II V CI we address 



insulating phases (Sec. IIV B|) 
more general questions like the behavior of one-particle 
correlations and scaling properties at finite temperatures. 

In harmonic traps we normalize using a length scale 
set by the combination lattice-confining potential, 



-1/2 



so that 



N 



^ j/=l 



(23) 



(24) 



In addition, instead of the density p = N b /N , relevant to 
the periodic or open case, we consider the characteristic 
density USUI 



N b a/(. 



(25) 



As shown in Ref. [36j up to p ~ 2.6-2.7 there is no Mott 
insulator in the trap. For larger values of p a Mott- 
insulating phase appears in the middle of the system. 



A. Superfluid case at T = 



noticed in n k around k = 0. For fc^T = O.lt [Figs, Elk) 
and Ef e )] , the energy is 3% larger than in the ground 
state and the peak in nk=o is less than one-third of its 
value at T = 0. For larger temperatures, like ksT = 0.5t 
in Figs. 0[c) andtTTf ) . almost no peak can be seen in n,k=o 
as compared with the one in the ground state. This is 
similar to the results obtained for the box in the previ- 
ous section, with a difference being that in the box the 
density distribution is not affected by the temperature. 




FIG. 7: Density (a)-(c) and normalized momentum distribu- 
tion function (d)-(f) of 200 HCB's in a trap with V 2 a 2 = 10~ 4 £ 
(p — 2) at different temperatures (solid line). The temper- 
ature [energy] of the system in each case is fcsT = O.Olt 
[E = -37.7H] (a),(d), O.lOi [-36.44i] (b),(e), and 0.504 
[—8.484] (c),(f). The dashed line in all the figures depicts 
the ground-state result (E = — 37.72t). 



In Fig. El we show density and momentum profiles in 
a trap with 200 HCB's (p = 2) for different tempera- 
tures (solid line) and compared to the ground-state case 
(dashed line). As for the fermionic in Figs.^d) Elf), 
the changes of the density profiles with the temperature 
in Figs. Et a )~EI c ) are the ones expected for fermions. 
[HCB's and fermions exhibit identical density profiles, 
Eq. (EI)-] For temperatures much smaller than the Fermi 
energy, which is of the order of t for these systems, the 
density profiles almost do not change. The same occurs 
with the total energy E of the trapped cloud, as seen 
from their values reported in the caption of Fig. El On 
the other hand, the behavior of nk, related to off-diagonal 
one-particle correlations, is very different to the one of 
the density. At k B T = O.Oli [Figs.El>) and Eld)], when 
the energy of the system has changed by 0.03% with re- 
spect to the ground-state energy, changes can be already 



Other quantities of relevance to the harmonically 
trapped case are the natural orbital occupations [Figs. 
Et a )"Et c )] an d the wave function of the lowest natural 
orbital [Figs. Etd)-|Htf )]■ In Figs. EtdHEff ) we normalize 
the natural orbital wave function following Ref. 

ip = R 1 ' 2 ^, R = (N b (/a) 1/2 . (26) 

Like rifc, the natural orbital occupations exhibit a very 
strong dependence on the temperature, which can be un- 
derstood since they are also related to the off-diagonal 
one-particle correlations. More interesting, and qualita- 
tively different to the case in the box, is the behavior 
displayed by the lowest-natural-orbital wave function in 
Fig.|SJ With increasing temperature, for large fillings, the 
weight of the lowest natural orbital in the middle of the 
trap decreases, and for fc^T = 0.5t [Fig.[Sff)] it is exactly 
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zero. This behavior of the wave function is accompanied 
by the appearance of a degeneracy in the occupation of 
the lowest natural orbitals. These two effects are very 
similar to the ones generated by the increase of the fill- 
ing in the ground state of the system and the formation of 
a Mott insulator in the middle of the trap j^j. However, 
as seen in Figs.Efa)-0c) no Mott insulator is created by 
an increase of the temperature. 
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FIG. 8: Natural orbital occupations (a)-(c) and normalized 
wave function of the lowest natural orbital (d)-(f) for the 
same systems shown in Fig. Q Nb — 200 and Via 2 = 10~ 4 t 
(p — 2). The temperature [energy] of the system in each 
case (solid line) is k B T = O.Oli [E = -37.7U] (a),(d), O.lOt 
[-36.44t] (b),(e), and0.50t [-8.48i] (c),(f). The dashed line in 
all the figures depicts the ground-state result (E = — 37.72t). 
The insets in (a)-(c) display in more detail the lowest 11 nat- 
ural orbital occupations. 



In order to understand the above effect it is important 
to realize that the spectrum of noninteracting particles in 
a combination lattice-harmonic potential |3|| |37], |33 |33 
is very different to the one of the harmonic oscillator in 
the continuum. (In the latter case one can intuitively 
realize that the maximum weight of a condensate, or of 
the largest eigenvalue of the one-particle density matrix, 
occurs in the middle of the trap.) In a lattice with a 
superposed harmonic oscillator, the eigenvalues of the 
noninteracting Hamiltonian reduce their weight in the 
middle of the system when their energy increases. For 
energies larger of 2t [for the Hamiltonian in Eq. ©] [Zfl ]. 
the eigenstates of the Hamiltonian start to be localized at 
the sides of the tr ap; i .e., they have zero weight in the cen- 
ter of the system |3(| |33, |3g, |3!j . At zero temperatures a 
Mott-insulating domain in the center of the trap signals 



that these states are populated [36j. At finite tempera- 
tures the occupation of localized states occurs even when 
there is no Mott insulator in the system, which explains 
why the lowest natural orbital can exhibit a behavior like 
the one seen in Figs.[3{e) and|Htf) m the absence of the 
insulating core. 

Before analyzing the finite-temperature one-particle 
correlations in the confined system, which explain the 
previously observed effects in and the natural orbital 
occupations, we present in what follows an example of 
the consequences of the temperature in a system that in 
its ground state exhibits a coexistence of superfluid and 
Mott-insulating phases. 



B. Mott insulator is present at T = 

In Fig. we show density and momentum profiles of 
a system with 300 HCB's (p — 3) for different tempera- 
tures (solid line) and compared to the ground-state case 
(dashed line) . As for the superfluid case discussed in the 
previous subsection, density profiles are almost not modi- 
fied for fcgT <C t. Increasing the temperature one can see 
in Fig. OUb) that as fcgT approaches t the Mott insulat- 
ing (n = 1) plateau in the middle of the trap disappears. 
The effects of the temperature in rik are also similar to 
the ones in the case with no Mott insulator, strongly 
depends on the temperature in the system. 




FIG. 9: Density (a),(b) and normalized momentum distribu- 
tion function (d),(d) of 300 HCB's in a trap with Via 2 = lQ~ 4 t 
(p — 3) at different temperatures (solid line). The tem- 
perature [energy] of the system in each case is ksT — O.lt 
[E = 158. 7t] (aj,(c) and 0M [185.9*] (b),(d). The dashed line 
in all the figures depicts the ground-state result (E = 157. At). 

It is worth noticing in Figs. Etc) and Eld) that even 
in the presence of a Mott-insulating phase, at zero tem- 
perature, nk exhibits a sharp peak at k = due to the 
superfluid phases at the sides 0, |23| ■ The effects of the 
Mott insulator in arc reflected by a large population 
of k states around ka = ±7r and an increase of the full 
width at half maximum of the k = peak. These are 
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characteristics of the system that remain at finite but 
very low temperatures. Increasing the temperature [Fig. 
Eld)] the peak at k = disappears. On the other hand, 
the high-momentum tails remain almost unmodified with 
respect to the ground-state case as they reflect the prop- 
erties of short-distance correlations, related to the den- 
sity profiles, which are much less sensitive to temperature 
effects. 
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FIG. 10: Natural orbital occupations (a),(b) and normalized 
wave function of the lowest natural orbital (c),(d) for the same 
systems shown in Fig. 03 N b = 300 and V 2 a 2 = 10~ 4 t (p = 3). 
The temperature [energy] of the system in each case is k B T = 
O.lt [E = 158.7t] (a),(c) and 0.5t [185.9t] (b),(d). The dashed 
line in all the figures depicts the ground-state result (T = 
and E = 157.4). The insets in (a),(b) display in more detail 
the lowest 11 natural orbital occupations. 



At zero temperature, the natural orbital occupations 
exhibit a clear signature of the presence of the Mott- 
insulating core in the trap. A plateau with A, ; = 1 is 
present, reflecting the existence of single occupied states. 
In this case the lowest natural orbital is degenerate due to 
the splitting of the system by the Mott-insulating core. 
Two identical quasicondensates can be observed at the 
sides of the Mott core [Figs. QUc) andHOTd)]. The in- 
crease of temperature reduces the occupation of the low- 
est natural orbital [they become more localized, Figs. 
HUT c) and ITUT dl] , but does not destroy their degeneracy. 
This degeneracy in absence of a Mott-insulating state is, 
as explained in the previous subsection, an effect that 
only appears at finite temperatures due to the popula- 
tion of localized states at the sides of the trap. Finally, 
one should notice that the plateau with = 1 disap- 
pears in Fig. HOf b) along with the disappearance of the 
Mott plateau in Fig.^b). 



C. Correlation functions and scalings 



are measured with respect to a fixed point Xj, while Xi 
is changed all over the system. In Fig. Illf a) the correla- 
tions are measured with respect to the middle of the trap 
and in Fig. II If b") with respect to two points at the sides 
of the Mott-insulating core present at T = 0. 




r / ' " \ ' 
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1-2-10123 

FIG. 11: One-particle correlations corresponding to the sys- 
tems analyzed in Sees. I1VAI and ITvPll (a) N b = 200, 
V 2 a 2 = 10" 4 i (p = 2), and x//( = 0. (b) N b = 300, 
V20 2 = 10~ 4 £ (p = 3), and Xj/C, = ±1.5. Different temper- 
atures are denoted by k B T — (dashed line), O.Olt (dotted 
line), O.Wt (solid line), and 0.50t (dash-dotted line). 

In the ground state, the one-particle density matrix 
de cay s as a power law pij ~ \x{ — x 3 -\~ 1 / 2 for < ni,rij < 
1 [23. The introduction of a small temperature = 
O.Olt can be already noticed in Fig. Illf a) as a faster 
decay of correlations at long distances. At temperatures 
larger than ksT = O.lt, for the system sizes of the figure, 
the one-particle density matrix decays exponentially to 
10~ 8 before reaching the borders of the trap. Due to 
the space varying density one can notice that in contrast 
to the box, in a harmonic trap one cannot see a single 
correlation length, which would mean a straight line in 
all the semilogarithmic plots of the figure. Still one can 
calculate the second moment of the one-particle density 
matrix £ [Eq. (|22|l ] as a sort of an averaged correlation 
length. 



TABLE I: Second moment of the one-particle density matrix 
(£), in harmonic traps (V2C1 2 — 10~ 4 t), for different temper- 
atures, and characteristic densities. We also present values 
of £ calculated in boxes (N = 500) with the same densities 
than in the center of the trapped case. Notice that in contrast 
to the characteristic density, the density in the center of the 
trap changes with increasing temperature. The two values of 
p correspond to ksT — O.lt and 0.5t, respectively. 



i/a 


(trap) k B T = O.lt 


k B T = 0.5t 




f/a (box) 


P 


= 0.5 


8.7 /9.7 


1.9 12.0 
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= 0.33,0.29 
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= 1.0 


10.4 /ll. 1 


2.1 /2.2 
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= 0.48,0.47 
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= 1.5 


10.5 110.3 


2.1 12.1 


P 


= 0.61,0.61 


P 


= 2.0 


9.7 /7.7 


2.0 11.6 


P 


= 0.75,0.75 


P 


= 2.5 


8.8 /2.8 


1.8 /l.l 


P 


= 0.93,0.87 


P 


= 3.0 


8.5 /0.0 


1.7 jO.7 


P 


= 1.00,0.95 



In Fig. 1111 we show the behavior of one-particle cor- 
relations with increasing temperature for the two cases 
analyzed in the previous subsections. Correlations (ptj) 



We present in Table ^ results for £ in harmonic traps 
for two temperatures and six values of p. To the right we 
show results obtained in boxes with densities chosen to 
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be identical to the ones at the center of the harmonically 
trapped cloud. One can see that the results in both cases 
are similar far from the region where the Mott insulator 
sets in the middle of the trap (p =0.5-2.0 in Table HJ, so 
that one can estimate £ in harmonic traps using results 
from a box. This is in agreement with recent results re- 
ported for other finite-temperature correlation lengths in 
trapped bosonic systems with no lattice |4l| . The reason 
for the agreement between £ in the trap and in the box is 
that in the first case £ is dominated by the contributions 
of the middle of the system, where the density has its 
maximum value and it is "relatively uniform." As one 
approaches n — 1 in the middle of the trap, or p = 2.6, 
the argument above fails (p = 2.5, 3.0 in TableQJ since the 
correlation length in the center of the cloud approaches 
zero (see £ vs p in Fig. when p — » 1) and regions with 
smaller densities start to dominate the value of £. For 
those cases an exact calculation of £, given the density 
profile, is required. 
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FIG. 12: Scaling of rik = o (thick lines) and the lowest natural 
orbital occupation (thin lines) vs Nb for a constant character- 
istic density p = 2. Different temperatures are plotted with 
the same convention of Fig. 1111 fcsT = (dashed line), O.Olt 
(dotted line), 0.05t (solid line), and O.lOt (dash-dotted line). 

The exponential decay of the one-particle density ma- 
trix implies that when the size of the system is larger than 
the averaged correlation length, the momentum distribu- 
tion function and the occupation of the lowest natural 
orbital stop changing with increasing system size. The 
size at which this occurs depends on the temperature, as 
the averaged correlation length decreases with increas- 
ing the temperature, and also depends on the charac- 
teristic density. In Fig. 1121 we show how rik = o and Ao 
scale at four different temperatures and starting from 
small system sizes (close to the ones achieved experimen- 
tally [20|1. At T = the increase of both quantities is 
~ s/Nhj reflecting quasi-long-range correlations present 
in pij |27| . For k B T = O.Olt, the departure from the 
zero-temperature values occurs when the trap has ~ 100 
HCB's. For k B T = 0.05i it occurs around N b = 20, and 
for k B T = O.lt even the smallest system with 10 HCB's 
is different to the ground state. 

To conclude this section we show in Fig. 1131 a compari- 
son between density and momentum profiles for 100 and 



400 HCB's at k B T = O.lt and p = 2. While even at zero 
temperature the density profiles do not differ [53, H(|, 
the ground-state peak nk=o would have been 2 times 
larger for 400 HCB's than for 100 HCB's. At k B T = O.lt 
both momentum distribution functions are almost indis- 
tinguishable. 




FIG. 13: (Color online) Density and momentum profiles in 
two traps with different filling, but the same p — 2 and fcsT = 
O.lt. 



V. GRAND-CANONICAL VS CANONICAL 
ENSEMBLE 

In the previous sections we have discussed effects of the 
temperature on trapped HCB's in ID. The starting point 
for our calculations was the grand-canonical ensemble, 
in which the system is assumed to be in thermal equi- 
librium with a large reservoir with temperature T and 
chemical potential p. The chemical potential was then 
chosen to obtain the desired average number of particles 
in the trap. In this section we analyze the changes intro- 
duced by the grand-canonical fluctuations of the particle 
number with respect to a fixed- Nb canonical description, 
which may be more relevant to describe trapped ultracold 
quantum gases where no particle reservoir is available. 

In the thermodynamic limit both descriptions are 
known to provide the same predictions |42j. On the 
other hand, for nonintcracting bosonic systems with a 
mesoscopic number of particles, it has been shown that 
the differences between the grand-canonical and canon- 
ical condensate fractions can be as large as 10% (for 
Nb = 100) close to the BEC transition point and decreas- 
ing logarithmically with increasing number of particles. 
In the present work we have been dealing with the oppo- 
site case — i.e., infinite repulsion. Interactions are known 
to suppress fluctuations of the number of particles in the 
grand-canonical ensemble |42^ . but since recent experi- 
ments with HCB's on optical lattices achieved only up to 
20 HCB's in around 50 lattice sites, it is useful to present 
an estimate of the difference between both ensembles for 
such small systems. 

In order to obtain the canonical one-particle density 
matrix we use the ground-state approach of Ref. |27j . 
We calculate the Green's function of all states \^ 7 hcb) 
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with Nb bosons in N lattice sites — i.e., of N s = N\/(N — 
Nb)lNbl states. The canonical Green's function at tem- 
perature T is obtained as the sum 

G % = zcT,^ En/kBT ^HCB\bS\n CB ), (27) 

71=1 

where e ~ E ^/ k BT - ls ^ e ener gy f state \^hcb)) * s ^e 
Boltzmann factor and Z c the canonical partition func- 
tion: 

Z c = J2 e ~ En/kBT - (28) 

n=l 

The canonical one-particle density matrix is then 
1 N ° 

n=l 

= Gg + (1 - 2Gg) . (29) 

In Fig. El we show results obtained for the grand- 
canonical (E) and canonical (E c ) energies of 10 HCB's 
in a box with 50 lattice sites as a function of the tem- 
perature. At the scale of the figure they are indistin- 
guishable. More information can be obtained in the in- 
set (a) where we plot as thin lines the energy difference 
between both ensembles as a function of the number of 
particles in boxes with densities N^/N = 0.2. As seen 
in this inset even for such small systems the difference 




FIG. 14: (Color online) Energy (in units of t) vs temperature 
(also in units of t) for a box with 50 lattice sites and 10 parti- 
cles. The grand-canonical (E) and canonical (E c ) results are 
plotted as solid and dotted lines, respectively. Inset (a) shows 
E-E c (thin lines) and SE = (E - E C )/\E C \ (thick lines) vs 
N b for k B T = 0.1t (solid line), 0.2t (dotted line), 0.5t (dashed 
line). Thin solid lines following SE exhibit 1/Nb behavior. In- 
set (b) depicts Sn^ for HCB's and 5nl for fermions vs Nb (see 
text). The temperatures are in this case fcsT = O.lt (solid 
line for Srik, dashed line for &n k ) and fcsT = Q.2t (dotted line 
for Srik, dash-dotted line for <5n{). Thin solid lines following 
these results exhibit I/JV5 behavior. In both insets the den- 
sity p = Nb/N — 0.2 was kept constant when changing the 
number of particles. 



almost does not change with Nb, and it is always smaller 
than the energy unit t. Considering that the modulus 
of the energy increases linearly with the system size, 
the relative difference between both ensembles decreases 
SE = (E — E C )/\E C \ ~ l/N b [thick lines in inset (a)]. 
For Nb — 10 and N — 50 one can see that SE is below 1% 
for temperatures up to ksT = 0.5t. As the temperature 
increases beyond fc^T = t the differences between E and 
E c start to decrease, which together with the decrease 
of the modulus of E c shown in Fig. El produces a sat- 
uration of SE at around 2% for fc^T > lOt. Then for 
Nb = 10 and N = 50 the maximum SE is just a 2% of 
the energy. We have also studied other densities keeping 
TV = 50, and the results obtained for the maximum SE 
were exactly the same 2%. 

While Kinoshita et al. [2l| used the energy of the sys- 
tem to confirm the achievement of the hard-core limit, 
Paredes et al. ppl ] considered the momentum distribu- 
tion function. In the inset (b) of Fig. El we show the rel- 
ative difference [Snk = (J2k \ n k~ n k D/Efc n k)} between 
the grand-canonical rik and canonical nj? calculation of 
the momentum distribution function. The relative differ- 
ences for nfc although larger than the corresponding ones 
for the energy are still small and also reduce ~ 1/Nb with 
increasing the system size. For Nb = 10 and N = 50 they 
are smaller than 1% up to ksT = 0.2t. 

It is also useful to calculate the differences between the 
grand-canonical and canonical ensemble for the equiva- 
lent noninteracting fermions. This may be relevant for 
systems like the ones recently achieved experimentally 
by Kohl et al. |44j . For noninteracting fermions, the 
energy differences between both ensembles are identi- 
cal to the ones of the HCB's due to the mapping, Eqs. 
(0) (J5J, so that as shown in Fig. E] and its inset (a) 
they are small. For the fermionic momentum distribu- 
tion function the HCB results do not apply. We have 
also calculated the fermionic relative difference [Snj. = 

Efc \ n k~ n k D/Efc n k )] between the grand-canonical 

f f c 

nt and canonical nt' calculation of n^. They are larger 

than for the ones of the HCB's as shown in the inset (b) of 
Fig. El However, they are still small for the experimen- 
tally accessible system sizes. For Nf = 10 and ./V = 50 
they are smaller than 3% for fc^T = 0.2i. Apart from an 
even-odd effect that decreases with increasing the tem- 
perature, Snj. also decreases ~ 1/Nb with increasing the 
system size. 

The introduction of a harmonic trap docs not (quali- 
tatively) change the results obtained in a box. In Fig. 1151 
we show the grand-canonical and canonical results of the 
energy in a harmonic trap with 10 particles and p = 2 as 
a function of the temperature. Contrary to the box, in 
a harmonic trap the energy is not bounded from above 
for very large temperatures. This is because the HCB 
cloud can increase its size and consequently its potential 
energy. The energy differences between both ensembles, 
when changing the number of particles keeping p = 2 
constant, are shown in the inset of Fig. El As for the 
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box they are almost independent of Nt,, and smaller than 
t. The results for the relative differences between the 
momentum distribution functions for HCB's and nonin- 
teracting fermions in both ensembles are also shown in 
the inset. Their behavior is very similar to the one of the 
box in inset (b) of Fig. 
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FIG. 15: (Color online) Energy vs temperature for a harmonic 
trap with 10 particles and p = 2. The grand-canonical (E) 
and canonical (E°) results are plotted as solid and dotted 
lines, respectively. The inset shows E — E° (thin lines) and 
Srik (thick lines) vs Nt for fcsT = 0.3/ (solid line) and 0.5t 
(dotted line). The <5n{ of the fermions is also shown as a 
function of Nj (thick lines) for ksT — 0.3/ (dashed line) 
and 0.5/ (dash-dotted line). Thin solid lines following the 
Silt and Snl results exhibit 1/Nt behavior. In the inset the 
characteristic density p — 2 was kept constant when changing 
the number of particles. 



As mentioned before, Herzog and Olshanii 43] dis- 
cussed the grand canonical and canonical differences be- 
tween the condensate fraction for noninteracting bosons 
in harmonic traps. At finite repulsive interactions, in 
ID, there is no BEC even at zero temperature. Still, for 
the HCB's we have calculated the differences between 
the largest eigenvalue of the one-particle density matrix 
(equivalent to the condensate occupation for BEC jU) 
in the grand-canonical and canonical ensembles. As for 
Silk in the inset of Fig. ED we find that the difference 
between them decreases ~ 1 /Nt with increasing number 
of particles in the system. This clearly contrasts with the 
~ l/ln(Ah) obtained for the noninteracting case 

We conclude by explicitly showing in Fig. 1161 the den- 
sity profiles and momentum distribution functions of 10 
HCB's in a harmonic trap with p = 2 at fc^T = 0.5t as 
obtained from the grand-canonical and canonical descrip- 
tions. They are basically indistinguishable. Then, even 
for the small system sizes achieved experimentally, one 
can rely on the grand-canonical description for strongly 
correlated HCB's for the physical quantities described 
here. The same conclusion applies to noninteracting 
fermions. 
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FIG. 16: (Color online) Density and momentum profiles of 10 
HCB's in a harmonic trap with p = 2 at fcsT = 0.5/. The 
grand-canonical and canonical results are plotted as solid and 
dashed lines respectively. 



VI. CONCLUSIONS 

We have presented an exact study of finite-temperature 
properties of HCB's confined on ID lattices. In or- 
der to solve this problem we have used the Jordan- 
Wigncr transformation to map HCB's into noninteract- 
ing fermions. After the mapping, properties of Slater de- 
terminants allowed us to obtain an exact expression for 
the HCB one-particle density matrix in terms of deter- 
minants of N x TV matrices, which are evaluated numer- 
ically. Our approach represents an alternative for finite 
systems to previous works that considered the thermo- 
dynamic limit |U 0, HE HE Hll f° r periodic and open 
chains and in which Toeplitz determinants were involved. 

We have shown that the effects of small finite tem- 
peratures are very important when dealing with quanti- 
ties related to off-diagonal one-particle correlations like 
the momentum distribution function and the natural or- 
bitals. These finite-temperature effects depend strongly 
on the system size. On the other hand, obscrvables re- 
lated to diagonal one-particle correlations (identical for 
fermions), like density profiles, are much less affected at 
low temperatures. Explicit results for the behavior of all 
these quantities versus temperature were given for system 
sizes that range from the ones recently achieved experi- 
mentally up to 20 times larger. 

Finally, we have compared grand-canonical and canon- 
ical results for energies and momentum distribution func- 
tions of HCB's and noninteracting fermions for small sys- 
tems, like the ones achieved experimentally. In spite of 
the mesoscopic number of particles we have shown that 
for these system sizes the effects of the grand-canonical 
fluctuations of the particle number are very small and 
one can rely on a grand-canonical approach. 

Although all our calculations are exact for infinite on- 
site U repulsion, for very strong but finite U our conclu- 
sions are still valid. In this case 1/U acts like a perturba- 
tion to the noninteracting spinless fermion Hamiltonian 
|45j . Rey et al. [^El have recently presented results ob- 
tained by exact diagonalization that support the above 
conclusion [39l | . A connection to experimentally relevant 
parameters can be also found in Rcf . j3^ . 
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